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ABSTRACT 


Recent developments in pseudo-holography provide the means to 
generate a three-dimensional image of an object from any radiation it 
scatters. One promising technique uses scattered radiation to cast 
shadows of a Fresnel zone plate, referred to as the coded image, 
which are recorded on film. The coded image is then decoded by 
passing coherent light through it which is diffracted producing an 
image of the radiation source. 

In this paper, a brief illustration of pseudo-holography using 
scattered thermal neutrons is presented, along with a development of 
a general expression for the coded image. The remainder of the paper 
is devoted to developing a scattered neutron flux source function by 
diffusion theory methods for a non-reentrant, steady state, homogeneous 
body of simple geometry. To test the accuracy of the source function, 
the computed scattered flux is compared with some measured values. 

The comparison supports the diffusion approximation for detection 


methods which do not have a strong forward detecting preference. 





I. INTRODUCTION 


Recent developments in pseudo-holography provide the means to 
generate a three-dimensional image of an object from any radiation it 
scatters. One promising pseudo-holographic technique is basically a 
two-step process. The scattered radiation is used to cast shadows of 
a Fresnel zone plate (hereafter termed the coded image) which are 
recorded on film. The coded image is then decoded by passing coherent 
light through it which is diffracted producing an image of the radiation 
Scattering source. This technique, first proposed for radiation sources 
by Mertz and ounce in x-ray astronomy, and extended to gamma ray 


ENC and Roeerece has been successfully applied to 


sources by Barrett 
actual image formation using x-rays scattered from an object by Sider, 
et sala" The use of scattered neutrons in image formation is presently 
under study by Rose, et ieee and provides the motivation for this paper. 
The pseudo-holographic technique using scattered neutrons differs 
from the use of scattered x-rays only in that the film and Fresnel zone 
plate respond to neutrons, not x-rays, in forming the coded image. It 
follows that an analytical representation of the coded image can be 
developed by describing the scattered neutron flux, modifying the 
scattered flux with a Fresnel zone plate response, and applying a film 
(detector) response function to the modified flux (shadows) seen by the 
film. The first step in this analysis is an accurate description of the 
uninterrupted scattered flux. It is the intent of this paper to develop 


an analytical method which describes scattered neutron fluxes emerging 


from objects in a beam of thermal neutrons. 





First, a brief illustration of pseudo-holography is presented 
along with a development of a general expression for the coded image. 
The remainder of the paper is devoted to developing a scattered flux 
source function by diffusion theory methods for a non-reentrant, 
steady state, homogeneous body of simple geometry. To test the 
accuracy of the source function the computed scattered flux is compared 


with some measured values. 





II. BACKGROUND 


A. Concept Illustration 


The basic concept of the pseudo-holographic technique noted above, 
as applied to thermal neutrons, is illustrated in Figure 1. A Fresnel 
zone plate is positioned relative to the scattering body and a suitable 
film-converter camera as shown. The Zone plate has a series of equi- 
area annular zones which are alternately transparent and opaque to 
thermal neutrons and consist, for example, of aluminum and cadmiun, 
respectively. Considering only a point scattering source of thermal 
neutrons, the display on the emulsion film (the coded image) is a 
shadow of the zone plate the shape and size of which depend monotoni- 
cally on the position of the source in the z-plane, (x,y); the position 


of the zone plate relative to the source plane, z,-z; and the position 


il 
of the camera relative to the zone plate, Zn—Z4+ When a beam of 
coherent light passes through the coded image, the light is diffracted 
into a focal spot which is an image of the original point scatterer. 
Considering any object to be a set of point scatterers, any specimen 
can be imaged. Each plane in the object will have a different focal 


plane when imaged, so that the three-dimensional information in the 


coded image can be viewed tomographically. 


B. Coded Image 


* 


Neglecting the detector resolution”®, assuming a linear film 


response, and noting the film output is different transparencies for 


*The detector resolution is idealized by considering the film and 
converter to be in the same plane and neglecting the film resolution. 
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points on the film, the coded image can be represented by the neutron 
flux shape in the x"y"-plane. As shown in Figure 1, (x",y') describes 
points in the film-converter plane at z = Zo 3 (x',y') describes points 
in the zone plate plane at z = Z43 and (x,y,z) describes points in 

the scattering body. 

Let Fi, y,2,2), the scattering rate density, represent the 
number of neutrons per unit time per unit volume per solid angle 
scattered in the volume dxdydz about x,y,z and in the direction dQ 
about 2. Let G(x,y,z,2 > x" yy" 42, 2") represent the angular flux 


Green's function to the film-converter plane, then 


-T (r, rE) 


GiGtey cer ae Se Weezer = = §/) _ 2. 6 (Q"=2) (1) 
9¥ 9S _dt »Y » q°s aie ae rr" ee 


where r,r' and r" are vectors from the origin (see Figure 1) to points 

in the scattering body (x,y,z), to points in the zone plate plane 

(x',y',z'), and to points in the film-converter plane (x",y"' 52). 

respectively; and n(r,r") is the "optical path length" between r and r". 
In this case it is convenient to consider n(r,r") in two parts, 


the “optical path length" in the scattering body, n, x"); and the 


"optical path length" outside the body, nN, (Eor")- Then 


_ TT = 1 = tt 
amr) 2 ony) | ony (er) (2) 
where the second exponential on the right hand side represents the Zone 
plate transparency which is defined by a set of step functions or a 


9 
series expansion . 





The neutron flux, B(x", y'sZy)s at the film-converter plane is 
then 


©) 


O(x",y",Z,)= | | | Fo (x,y,2,2)6(x,y,2z,2 > NY's Z, G2") dao" dx 
2 Q"Vol 


If an angular emitted neutron current from the surface of the body can 


be described as K(x, y,z,%), then 


(4) 


Gay" »Z,)= 


I 
IQ-—. 


| K(xsy,2 G(x, y,2z,& = x" sy" »Z, 2°) dQdQ" dr 
S 


where | represents integration over the body surface, dr. is a surface 
ee oe the body and in this case n, (x>r") = 0. Note that with 
the zone plate removed, nN, (45°) = Q in equations (3) and (4) above, 
and the uninterrupted scattered neutron flux at (x",y") in the film- 
converter plane is described. 

To evaluate equations (3) or (4) an expression for the emitted 
neutron current or the scattered neutron source must be developed via 
a suitable neutron transport approximation. The accuracy of the source 
function can be verified by computing the uninterrupted scattered flux 
from equations (3) or (4) and comparing its shape with measured values. 
As noted above, the flux shape in the (x",y") plane is sufficient to 


represent the coded image. 
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IIIT. SCATTERED FLUX FROM A REGULAR SOLID 


For the development of a coded image of the interior structure to 
be practicable, the thermal neutron scattering cross section, a of the 
scattering object must be much greater than the thermal neutron absorp- 
tion cross section, ae This suggests the possible application of the 
diffusion approximation to neutron transport in developing a scattered 
source function as used in equations (3) or (4), where some approximation 
will be required to define the directional dependence of the scattered 


source function. 


A. Scattered Source Function 

Breas assumed that an adequate test case is to examine the 
scattered flux from a regular solid in a beam of thermal neutrons 
of intensity I, (neutrons per unit time per unit area) as shown in 
Figure 2, where the composition of the solid is homogeneous and such 
that Z_>> De As mentioned earlier, either a scattering source function 
or an angular emitted current source function can be used in computing 
the scattered flux. To develop a scattering source function, the flux 
in the solid, (x,y,z), is derived by diffusion theory methods, and 


the solid is considered an isotropic volumetric source of strength 


uO (%, 952) 


An . The scattered flux at some point away from the solid, 
(X95 22) in Figure 2 would be the uncollided contribution from each 
point in the solid. An angular emitted current source function on the 
surface of the solid is defined directly from the emitted current 


derived via diffusion theory, where some emission angular dependence 
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must be assumed. The use of diffusion theory in computing reflection 
coefficients and the relative ease in computing the scattered flux from 
an angular emitted current Source suggests the use of the angular 
emitted source. 


The diffusion equation to be solved is 
DV ®(x,y,2) - B,O(x,y,z) + S(x,y,z) = 0 (5) 


where D is the diffusion coefficient of the solid and S(x,y,z) is an 

isotropic volumetric source. It is convenient, in terms of describing 

suitable boundary conditions, to let S(x,y,z) represent the first 

collision scattering density, i.e., 

S(O nee ae (6) 

So 

where y, is the total thermal neutron cross section of the solid and it 

is assumed that the first collision scattering density is isotropic. 

The boundary conditions are that ®(x,y,z) is zero at the extrapolated 

boundaries of the solid. 

The solution to equation (5) (see Appendix A) obtained by the 


method of eigenfunctions for the geometry in Figure 2 is 


11x ; miry ; nz 
O(x,y,z) = ) A sin 235] sin = sin ze (7) 
ea 1lmn c a b 
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where 


a2 S Taek 
S Oo 


i =o 22 Se an 
Im “ E (4L°By)  Blmm E A Fe 


_ -_b = _,_ |[mn(b-2D) nt _ nT (b—2D) 
K=e it - an sin [ea tbe20) aa * cos [a=20>]| 


= ; nim2D ni apse | 
+ d. sin [ 5 ] 1 a cos (pee 
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Noting the symmetry in Figure 2, x = c/2 and y = a/2 is the center 
of symmetry in the xy-plane, and considering points (X 0 9¥ 22) which do 
not see the ends of the solid, i.e., 2D<z_<b+2D, then the emitted 
current need only be defined for faces 1 and 2 of the solid. The emitted 


current from face l, i (c+2D,y,z), by Fick's Law is 


“- Zp : 
i (orao.y.2) = Meetpae 2 | attnad | ie 
x = ct+2D 
where 
d0(x,y,2) E: y 1T A aoe 1m (c+2D) Petey se 
ox Cc imn Cc 
x = ct2D I1mn 





niZ 
sin kz: ] 





bal 


and the emitted current from face 2 is 


x Gap iy ye ee a B 9®(x,y,z) 


dy (9) 
y = at2ZD 
where 
dO(x,y,2) _ mT mm (at+2D) 11x 
aa = > ee AL an COS ee | Sin [ar | 
y = at2D a 
sin [S| 


B. The Scattered Flux 

To use equation (4), with n(r,r") = 0, in calculating the 
uninterrupted scattered flux, O(x 9¥ 522) at (XY 22s it is necessary 
to define some angular dependence to the emitted current. Realizing 
it is difficult to make any generalization about the angular dependence 
of the emitted current, two possibilities for defining an angular 
dependence in terms of isotropy are: 

1. Replace the emitting face with a plane semi-isotropic source 
(hereafter termed the semi-isotropic emitted current assumption) of 
strength NGaaee) defined by equation (8) or (9). 

2. Assume the current density, j(x,y,2,2), is semi-isotropic 
(hereafter termed the semi-isotropic current density assumption) 


then 


ce (x,y,2) = [ a. jx,y,2,2)d2 (10) 
27 





2 


where j(x,y,z,%) = CQ and n is the unit vector normal to the emitting 
area. 
Based on the semi-isotropic emitted current assumption the angular 


emitted current from face 1, KK. 6ct2D,y,z,%) » is 


J” (c#2D,y,2) 


K_,(c+2D,y,z,2) = oT (11) 


sl 


and the angular emitted current from face 2, K 7 (%,at2D, 2%) , is 


7 Cen 
ee eee 
K 9 (x,at2D,z,&2) a (12) 
Based on the semi-isotropic current density assumption, the 
angular emitted current from face 1, K. 9 (c+2D,y,2,%), is 
K_, (e+2D,y 2,2) =n. j(ct2D,y,z,) = Cy cos ¢ C13) 


where cos ¢ =n. & and ¢ is then the angle between the neutron direction 
and unit normal to the emitting face. Evaluation of equation (10) 
yields 


J” (c+2D,y,2) 


oa T (14) 


Similarly the angular emitted current from face 2, Ki (x,at2D,z,%), is 


J" (x,a+2D,2) 
=» iia eee 


K 2 (%,at2D,z,2) = = cos > (15) 
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Using equation (4), with n(r,r") = 0, the uninterrupted flux, 
(x s¥522,) at (x s¥ 92,5 recalling that 2D<z_<b+2D, is determined 
for the two angular emitted currents as follows: 
1. For the semi-isotropic emitted current assumption: 
Age pet 2D<y <atzD and x 7c+2D then only face 1 contributes 
to the flux at (x s¥592,) and evaluation of equation (4) over 


Q and 2" yields 


J. (e+2D,y,z)dydz 
O(x VY _ 92 ) = as 
Oo 0 0 


Lae oe oe (16) 
2D 2D 2 [[,~ (e+20)] +(y -y)"+(2_-2) | 
bow cimillarly, if 2D<x_ <c+2D and y pared then only face 2 
contributes to the flux at (x0 s¥ 5920) and 
ctZD b+ZD ais 
5 ‘a Jj (x,atZD,z)dydz 
(X 0 5¥ 59%, a Sun PRINS Sara EEEE NOUN (17) 
a 2p 2h (x0 +[y -(at2D)] +(z -2) | 
Ce er y,7at2D and x7 ct2D, then faces 1 and 2 contribute 
to the flux at (xX s¥ 092) and 
O(x V5 2%) = equation (17) + equation (16) (18) 


2. For the semi-isotropic current density assumption: 
a. If only face 1 contributes to the flux at (XY 092,) the 


evaluation of equations (4) over 2 and %" yields 
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at+2D b+2D mt 
Mee ah ee cos >, J, (c+2D,y,2) dydz es 
Oo Oro 


2D 2D n| [x,~(ct20)] “Hyg-y) He -2)” | 


where the integration of equation (4) over &% requires that 


a x ~(c+2D) 
cos $, = cos |e py| = a ae ae ee 
, ae | Fx, (e#20) saa) +(z-z) lee 


b. If only face 2 contributes to the flux at (XV 9%) then 


c+2D b+2D a 
O(x ,y_,z ) = cos $,J (x,at2D,z)dydz 
ene 2 2 7) (20) 
a ee iD (x,-» +[y -(at2D)] “+(z_-2) | 
where 
y —(a+t2D) 
O 
cos %, = 


2 2 21/2 
I(x, x) +ly,, (a+2D)| +(z Zz) | 
c. Lf faces 1 and 2 contribute to the fiux at (Xo 5¥0 92) then 


(x99 59%) = equation (19) + equation (20) (21) 


To evaluate the above integrals, a numerical technique is 
conveniently applied. The series in equations (8) and (9) are uniformly 
convergent so that the series can be integrated term by term. Thus each 


term in equations (16), (17), (19) and (20) can be integrated numerically 





15 


and a suitable convergence established. Results of scattered fluxes 
for a specific regular solid computed as outlined above are contained 


in the next section. 
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IV. METHOD OF MEASUREMENT AND CALCULATION 


A. Method of Measurement 

Measurements of the scattered flux were obtained with an arrangement 
simulating the geometry of Figure 2. As shown in Figure 3, a regular 
solid specimen, a polyethylene (C,H, ) block (2.86 cm square x 5.16 cm 
long), was placed in a beam of thermal neutrons. A BF detector was 
used to count scattered neutrons. The detector longitudinal axis 
remained perpendicular to face 2 and in the yz-plane at x = 1.43 cm + 
2D, where D is the diffusion coefficient in cm. The center of the 
detector window was located at the following points (XY 0920) for 
measurement: 

1. (1.43 cm + 2D, 7.86 cm + 2D,z) where Zz was incremented by 
By.) cm from .5 cm + 2D to 4.5 cm + 2D. 

2. Same as 1, except eg iZ.o0 cm + 2D. 

The Pennsylvania State University TRIGA Reactor (Breazeale 
Nuclear Reactor Facility) provided the source neutrons. A D0 tank 
between the beam port entrance and the reactor core provided sufficient 
thermalization such that the thermal flux at the beam port exit at a 
reactor power of 1LOOKW was on the order of to! neutrons per one per 
second. Collimation was achieved by placing in the beam port two iron 
plugs with 1-1/2 inch holes followed by a copper plug with a 3/4 inch 
hole. Smaller collimations would reduce background scattering, however, 
to maintain the specimen totally within the beam the above sizes were 


the minimum allowable. The specimen was placed 15 inches from the beam 


port exit and at this point, it was estimated that the beam had diverged 
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to approximately twice the exit diameter with a flux on the order of 
10° neutrons per ant per second. The cadmium ratio in the beam was 
approximately 50. 

The BF, detector was 14.5 inches long with a 1 inch outer diameter 
and was filled to a pressure of 60 cm Hg with BF, 96% enriched in ae 
The detector was shielded with a cadmium jacket, and the window had a 
cadmium cover with a 1/16 inch diameter hole. By placing a neutron 
shutter between the beam port exit and the detector, an average signal 
to noise ratio of 2 was achieved. As shown in Figure 3 the neutron 
counting arrangement was typical for a gas filled proportional counter. 

The polyethylene specimen was held in place by an aluminum clamp 
attached to a horizontal bar 1/4 inches in diameter. Three of the 
side faces were covered with a cadmium sheet which reduced background 
scatter from the aluminum holding stand, and maintained the integrity 
of the non-reentrant boundaries. 

One minute background counts were taken at each point described 
above with only the specimen removed. The gross counts were also 
taken for one minute and were sufficiently high to yield counting 
fractional standard deviations on the order of .025. The counts 
obtained for points 5 cm Cs 7.86 cm + 2D) and 10 cm Ora 12.86 cm + 
2D) from face 2 were normalized with the z5 7 2D = 2.5 cm point. The 
results are displayed on Figures 4 and 5. As discussed earlier, the 
normalized curve represents the relative scattered flux shape as 
amended by the detector response. 

Additional measurements were made at 2.5 cm oes 5.36 cm + 2D) 


from face 2. The measurements were repeated with various strips of 


cadmium on face 2 as follows: 
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1. A 1 cm wide strip beginning at x = 1.5 cm + 2D and ending at 


x = 2.5 cm t+ 2D. 
2. A 2 cm wide strip beginning at x = 1.5 cm + 2D and ending at 
x = 3.5 cm + 2D. 


All counts obtained were normalized with count obtained at (1.43 cm 
+ 2D, 5.36 cm + 2D, 2.5 cm + 2D). The results are displayed on Figures 


oe sand &, 


B. Method of Calculation 

The method of calculation was outlined in the previous section, 
and in this case equations (17) and (20) applied. To compare the 
calculations with the measured values, a response function was inserted 


representing the BF, detector response. The response function was 


3 


determined as follows. 
Neglecting the finite area of the detector face, the probability 
of detection for a neutron entering the detector at an angle 9 with 


the detector face is 


1 - exp(-E r/cos 8) (22) 


10 


= -B 
where Le = Le is the thermal neutron absorption cross section of B 


at 60 cm Hg, r is the radius of the detector and r/cos 0 represents 
the path length of the neutron in the detector. For the geometry in 


Figure 3 


172 
(x ene ae 
cos § = - = (23) 


(x, -x) “+ ly -(a+2D)| “+(2-2)° 
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Considering the finite area of the detector window, a unit flux 
in the @ direction to the detector face produces a current of sin 6 


entering the detector, where from equation (23) 


1/2 
[yoo (at+2D)] = 


2 2 2 (24) 
(x -x) +[y,-(at2D) | +(z -z) 


sin 9 = 


The detector response function, R(X, oY 59% 2X9 22) » is then the fraction 
of the flux seen, equation (24), times the probability of detection, 
equation (22). 

Assuming a semi-isotropic emitted current and accounting for the 


detector response, equation (17) was modified to 


c+2D b+2D 


=" 
= R(X _sY¥_ 52 sXs¥5Z) J (x,at+2D,2z)dxdz 
O(X 9¥ 5225) = ondels ale (25) 


2D 2De2n | x.) “+ [y,- (at 20) +(2,-2)° ] 


Assuming a semi-isotropic current density equation (20) was modified to 


c+2D b+2D (26) 


++ 


2D 2D 1 (x,-0 “+ [y,-Cat20)] “+(25-2)"] 


Equations (25) and (26) were solved, as noted previously, by integrating 
each term numerically using as IBM system 370 Mod 168 computer and 
10 


Argonne National Laboratory subroutine DROMB ~. The series was truncated 


when the higher mode integral terms were less than .1% of the existing 





26 


sum. The thermal neutron diffusion properties of polyethylene were 
obtained from the ant-58007t, Calculations were made for each point 
measured by the detector and were similarly normalized. The results 
for the semi-isotropic emitted current and for the semi-isotropic 
current density are displayed on Figures 4 and 5. 

The calculations were repeated for the detector at 2.5 cm from 
face 2 as used in obtaining the measured values. Calculations were 
made at each point measured by the detector and were similarly 
normalized. Calculations with the cadmium strips on face 2 were made 


the semi-isotropic current density assumption only. The results are 


displayed on Figures 6, 7, and 8. 
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V. ANALYSIS OF RESULTS 


In analyzing the results on Figures 4 through 8 it is helpful to 
bear in mind that: 

1. The detector had a forward counting preference which decreased 
as the detector was moved away from the emitting surface (see equations 
(22) and (24)). 

2. The emitted current, used in determining the angular emitted 
current, resulted from a diffusion calculation. 

3. The angular emitted current was determined from assumptions 
regarding angular dependence. 

As shown by Figure 4, at 10 cm from the emitting face the measured 
flux shape was generally well represented by the calculated flux shapes 
with differences only at points near the boundaries of the emitting 
face. Coupled with the detector forward counting preference, the 
boundary errors would be expected from a diffusion calculation. The 
differences, as displayed in Figure 4, indicate that the boundary 
errors were the result of calculated emitted currents which were low 
near the edges of the emitting face. 

Additionally, a portion of the observed differences can be 
attributed to the assumptions used in defining an angular emitted 
current as the calculated flux shape for the semi-isotropic current 
emission differs from the calculated flux shape for the semi-isotropic 
current density emission. To quantitate somewhat, the maximum 
difference between the calculated relative flux for a semi-isotropic 
emitted current and the measured relative flux occurred at Presi 2D = 


4.5 cm and was 15.1%. The difference between the calculated relative 
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flux for a semi-isotropic current density and the measured relative 
flux at the same point was 18.3%. Thus, the errors displayed by 

Figure 4 were due primarily to errors in the calculated emitted current 
near the boundaries. 

As shown by Figure 5, when the detector was 5 cm from the emitting 
face, the boundary errors increased. This would be expected as the 
detector preference for points directly in front increased and the 
calculated emitted current was in error at the boundaries. Moreover, 
errors in the assumed angular dependence of the emitted current 
increased. To quantitate again at ae 2D = 4.5 cm, the difference 
between the calculated relative flux for a semi-isotropic emitted 
current and the measured relative flux was 26.8%. The difference 
between the calculated relative flux for a semi~isotropic current 
density and the measured relative flux was 33.82. 

As shown by Figure 6, when the detector was 2.5 cm from the 
emitting face, the differences increased further as would be expected. 
The differences in the calculated values for the semi~isotropic current 
density and the semi-isotropic current were small compared to the total 
errors. The principal source of the errors was again the calculated 
emitted current. Moreover, it should be noted that the block dimensions 
were only 1L x 2L, where L represents the diffusion length. 

The measured and calculated flux shapes at 2.5 cm from the emitting 
face with 1 cm and 2 cm cadmium strips on the emitting face are 
displayed on Figures 7 and 8, respectively. The boundary errors dis- 
played by Figure 6 are not as evident here, due to the normalization. 


In this case a significant portion of the face which had a reasonably 
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correct calculated emitted current was eliminated by the cadmium strip. 
As a result, the calculated flux at the Zz, 7 2D = 2.5 cm point had a 
larger contribution from the boundary areas yielding a low value. 

When normalized to this low value, the low values near the edges of 
the emitting face were dampened, decreasing the errors in the 
calculated relative flux. 

The above comments can be summarized as follows. Differences 
between the calculated values and the measured values were due primarily 
to boundary errors in the calculated emitted current, with very 
minor effects due to the arbitrary assumptions used in defining the 
angular dependence to the emitted current. As the detector forward 
preference was decreased, the observed differences decreased and the 
calculated flux shapes represented well the measured values. Although 
large differences existed between the calculated and measured values 
when the detector had a strong forward preference, measured dips in the 
relative flux over a strong absorbing strip were represented by the 


calculated values. 
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VI. CONCLUSIONS 


A. Summary 


The diffusion approximation was applied in calculating emitted 
neutron currents from a homogeneous regular solid in a beam of thermal 
neutrons. The calculated emitted current, with an assumed angular 
dependence, was used in computing the scattered flux shape about the 
solid. Calculated results, modified for the response of a cylindrical 
gas filled proportional counter, were compared to values measured with 
a BF, proportional counter at points which received flux contributions 
from only one face of the solid. The observed differences in the 
calculated and measured flux shapes were due primarily to errors in 
the calculated emitted current about the edges of the face, with very 
minor effects due to the assumed angular dependence of the emitted 
current. As the directional preference of the detector decreased, 
the observed differences decreased. Moreover, regardless of a signifi- 
cant directional detector preference, trends in flux shapes over strong 
absorbing strips were well represented. It can be concluded that for a 
detector with a lesser directional preference, such as a film-converter 
detector, the calculated flux shape would better represent the 
measured flux and would be sufficient to complete the analytical 


representation. 


B. Suggestions for Future Work 


The description of the scattered neutron flux shape represents 


the first step in an analytical representation of the pseudo-hologram. 
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Future work is suggested to complete the analytical representation by 
including the response functions for the Fresnel zone plate previously 
described and for a typical film-converter detector. The detector 
response function will have to be suitably approximated as was done 
herein for the BF, detector response function. The flux shapes so 
determined could then be compared to calculated flux shapes based on 
an ideal detector and an ideal Fresnel zone plate in an effort to 
develop a correction factor for pseudo-holograms practically obtained. 
One possible application of a correction factor would be in correcting 
pseudo-holograms which are decoded numerically on a computer ayeCenuae 
The correction factor, depending on a general application, would at 


least decrease some of the resolution loss experienced in the real 


system. 
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APPENDIX 


SOLUTION OF DIFFUSION EQUATION 


The geometry is as shown in Figure 2. Then the extrapolated 


boundaries 
c 
a 
b 


where c, a 


are 


ec + 4D 


a+ 4D 


bea 


and b are the block dimensions. 


The equation to be solved is 


DV8(x,y,2) - E.O(x,y,2) + S(x,y¥,2) = 0 (A-1) 


where the source is 


= - 5 (z-2D) 
Danes for 2D<z<b-2D (A-2) 


S(x,y,2) = 


and the boundary conditions are (x,y,z) extrapolated _ 


0 otherwise 


boundary : Se 


To solve (A-1) use 


Vo(x.y.2) + BL O(x,y5z) = 0 (A~4) 


to generate a linear independent, complete, orthogonal set of functions 


£5 6% 992) } for l,m,n = 1,2,3,++*+ where 
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bs an = sin (= *)si in) sin ne) and (A-5) 
2 2 Z 

Zee eT mT nT 

Lee = Ce) a {s ) + (F ) (A-6) 

Then 

S(x,y,Z) = ) 34am Pin X29 2) and (A-7) 
lmn 

(x,y,z) =) AL oO, (x,y,z) (A-8) 
1mn 


From the orthogonality property of to, (2722) and equations 


(A-2) and (A-7), then 


320 an 
S = aa a 
im = §1mn “(I +@ S| Oe) 


where 


eo gvetl [-:,. = ree) - Bl cos pe j . 


n7mZ2D 


nn2ZD ) 
b (A-10) 


Le sin(-—=— : ) eee a ” cos ( 





Substituting equations (A-7) and (A-8) into equation (A-1), using 


equation (A-4) and the linear independence of {$, (eV eZ} then 


a 
2 


B (itl 


A 
Imm” By an) (A-11) 
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where 
2 -_ 
L = D/L, 


Then using equations (A-8), (A-9), (A-10) and (A-11) 


VCRs) = 2 A 95 an $%2 92) (A-12) 
mn 


lmn 
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